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^O . We study gravitational collapse in higher dimensional quasi-spherical Szekeres space-time for 

matter with anisotropic pressure. Both local and global visibility of central curvature singularity 

has been studied and it is found that with proper choice of initial data it is possible to show the 

validity of CCC for six and higher dimensions. Also the role of pressure in the collapsing process 

, ■ has been discussed. 
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(N ; I. INTRODUCTION 

Usually, for cosmological phenomena over galactic scale or in the smaller scale, it is reasonable to consider 
f*^ ' inhomogeneous solutions to Einstein equations. Szekeres' [1] in 1975, gave a class of inhomogeneous solutions 
Q\ ' representing irrotational dust. The space-time represented by these solutions has no killing vectors and it has 

invariant family of spherical hypersurfaces. Hence this space-time is referred as quasi-spherical space-time. 

Recently, Chakraborty et al [2] have extended the Szekeres solution to (n + 2) dimensional space-time and 
^2 . generalized it for matter containing heat flux [3]. 

In classical general relativity, one of the challenging issues is gravitational collapse. This problem became 

jj important after the formulation of famous singularity theorems [4] and Cosmic Censorship Conjecture (CCC) [5]. 

Also in the perspective of black hole physics and its astrophysical implications, the end state of collapse (black 

5h . hole or naked singularity) is interesting. As there exists no formal method to address this problem so it is natural 

to study various examples of collapsing system (namely, Tolman-Bondi-Lemaitre (TBL) spherically symmetric 

£> model [6-15] or quasi-spherical Szekeres' model [16-18]) with a view to gain some insight. In general, these 

studies conclude that the local or global visibility of the central curvature singularity depends on the initial data. 

The above studies are mostly confined to the dust model, there are very few works on anisotropic stress 
[19] (confined to TBL model). For the last few years, there are attempts to study collapse dynamics in TBL 
model with anisotropic pressure to address the question "can non-zero pressures within a collapsing matter 
cloud avoid a naked singularity forming as the end state of a continual gravitational collapse?" But so far, the 
actual role the pressures play in determining the end state of collapse is not yet clearly understood. Due to 
complicated nature of the space-time geometry there are no works on gravitational collapse with anisotropic 
stresses in quasi-spherical model except recently by Chakraborty et al [20] where they have shown the role of 
pressure in 4 dimension. 

In this work, we extend this study to [n + 2) dimensional Szekeres' model and examine the role of the 
dimension on collapse dynamics. The paper is organized as follows: Higher dimensional Szekeres' model is 
described in section II, while collapse dynamics including the study of geodesic is presented in section III. In 
section IV, there are discussion and concluding remarks. Finally at the end there are two appendix dealing 
with detailed calculations. 
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II. HIGHER DIMENSIONAL SZEKERES' MODEL 

The metric ansatz for (n+2)dimcnsional Szekeres' space-time is of the form 



ds 2 = dt 2 - e 2a dr 2 - e 2f3 £ dx 2 (1) 

i=l 

where the metric cocfhcients a and f3 are functions of all space-time co-ordinates i.e., 

a = a(t,r,xi, ....,!„), (i = /3(t,r,xi, ....,x n ). 
Now Considering both radial and transverse stresses the energy momentum tensor has the structure 

T^l = diag(p, -p r , -p T ,-p T ) 
and the compact form of the Einstein equations are 

P=£r 

Pr = "^ (2) 

P T =Pr + T%T 



where F{r, t) = r jR n - 1 e ( - n+1 >(R 2 - f(r)) and C = e 



Further, the expressions for the metric functions are 

e (i =R(t,r) e^.zi >-,*«) j 



(3) 
e a = R> + Rv' 



and the evolution equation for R gives 



RR H — (n — 1)R 2 H — -R 2 = f{ r )i (/( r ) = arbitrary separation function) (4) 

2 n 2 



Also the function v satisfies 



e- 2 " 



2 ^[(n-2)^+2, IjIj ]=n(/(r)-l) (5) 

which has a solution of the form 

n n 

e-^AM^^ + ^B^rK + CW (6) 

with the restriction, 

n 

Y J B 2 -AAC = f(r)-l (7) 

i=i 

for the arbitrary functions A(r), Sj(r), (i = 1,2, ..,n) and C(r). 



As we are considering quasi-spherical gravitational collapse so it is natural to assume the initial configuration 
(from which the collapse has started) to be smooth everywhere. Thus p r should be regular initially at the center 
and blows up at the singularity. So a natural choice for p r is 



Pr = 



9(r) 
R l 



(8) 



where the arbitrary function g(r) has the form r l near r = to make initial p r finite (non-zero) at the centre 
r = and I is any constant. Hence, the expressions for matter density and tangential stress become 



P = 



H' + {n + \)Ev' 
R n (R' + Rv') 



and 



9(r) 
R l 



1- 



IR' 



n(R' + Rv 1 ) 



9'(r) 



nR l ~ 1 {R' + Rv') 



(9) 



(10) 



where H(R,t) = D{r) — /^^T+i) -^™ ' +1 ' (^ 7^ ( n + 1)) an d D(r) , an arbitrary integration function. 

Now, due to this choice of p r (see eq (8)) the evolution eq (4) for R can be integrated once and the radial 
velocity of collapsing shell at a distance r from the centre is given by 



R 2 = f(r) 



2H(R,t) 

nR 



n-l 



(11) 



This is termed as the equation of the collapsing process. 



III. COLLAPSE DYNAMICS 

To characterize the nature of the singularity (black hole or naked singularity), the event horizon of observers 
at infinity plays an important role. But formation of event horizon depends greatly on the computation of null 
geodesies whose computation are almost impracticable for the present space-time geometry. So a closely related 
concept of a trapped surface (a space-like 2-surface whose normals on both sides are future pointing converging 
null geodesic families) will be considered. Thus, if the 2-surface S r< t (r =constant, t ^constant) is a trapped 
surface then it and its entire future development lie behind the event horizon provided the density falls off fast 
enough at infinity. Hence mathematically, if K ^ denotes the tangent vector field to the null geodesies which is 
normal to S r t then we have 



K^ K» = 0, K» v K v 







Now the null geodesies will converge (or diverge) if the invariant K^. < (or K^. > 0) on the surface 

S r .t=0- 



As a consequence, it is easy to show that the inward geodesies converges initially and throughout the collapsing 
process while the outward geodesies diverges initially but becomes convergent after a time t a h{r) (time of 
formation of apparent horizon) given by 



R(t ah (r),r) = -y/l + f(r) 
Now using equations (8) and (11) we have 



g(r)R n+1 - [ (t ah (r),r) + -{n + 1 - /)i?"- 1 (t aft (r), r) - (n + 1 - l)D{r) = 



(12) 



From Appendix II, it is to be noted that the central singularity (at r = 0) forms at time to while a trapped 
surface is formed at a distance r at time t a h(r) and their difference is given by the equation (37). Thus if the 
trapped surface is formed at a later instant than to then it is possible for light signals from the singularity to 
reach a distant observer. Hence, t a h{r) > to is the necessary condition for formation of naked singularity and 
on the othcrhand, t a h{r) < to is the sufficient condition for black hole formation. Also it should be mentioned 
that this criterion for naked singularity is purely local. 

Further, as the time difference equation (37) is complicated, so to make a comparative study between t a h and 
to one can choose for simplicity I = (n + l)/2 and equation (37) takes the form 



tahir) - t = 



2n D g 1 - D 1 g - gw/DoJ D 



n+l 



fill 



r + 0{r 2 ) 



(n + lXgoVA^A) - £-g (V^o" + yjD - ^So) 
The following table shows the possibility of naked singularity or a black hole under different conditions 



(13) 



TABLET 



Choice of the parameters 


Naked Singularity 


Black hole 


(i) 3 i>0,£i<0 


Always possible 


Not possible 


(ii) 9l < 0, Di > 


Not possible 


Always possible 


(iii) 5i > 0, D x > 






737>^(i + V / i-4t^) 


«r<^( 1 + y /l -4r^) 


(iV) 9l < 0, Di < 






|^|<B±i(l + A / 1 __^^) 


m> i ¥( i +\/ i -£i%) 



From the table, to make the initial density gradient to be negative at the centre (i.e., p\ < 0) one must have 
{D\ — —+i) < (for i/_j > —I). For the first case (i.e., g\ > 0,Z?i < 0) p\ is negative definite and there is 
always naked singularity as in the dust model. Similarly, p\ is positive definite in the second case which leads 
to black hole solution same as dust model. For the third and fourth cases (when g\ and D\ have same sign) 
either NS or BH is possible depending on the restrictions given in the table I (see also figs 1-6). However, in 
the last two cases, for p\ > only black hole solution is possible but for pi < 0, both NS and BH are possible. 



Further, if we assume that D\ = = g\ then the time difference in eq (37) becomes 

n — 1 1 

l Dr 1 



tah{r) -t = x j- 



3re-l n+l 1 I IDn 

2 \ n 



y o 



3ti + 3-2; 

r ,,-1 



(14) 



We see that if n > 3 (with I < n + 1) then the first term on the right side will be dominating compare to 
other terms and hence we always have t a h < to- Thus in this case black hole is the only final state collapse 
for six and higher dimensional space-times. This distinctive result is similar to dust collapse [13, 16] and we 
conclude that CCC is valid in this case for six and higher dimension with anisotropic pressure. 



A. Study of Geodesies 



In this section, the nature of the singularity (NS or BH) is examined by studying the geodesies from the 
singularity. In particular, it will be investigated whether there exist one or more radial outgoing null geodesies 
which terminate in the past at the central singularity. For simplicity of calculation only marginally bound case 




Fig.l 



Fig. 2 





Fig.3 



Fig.4 





Fig. 5 



Fig.6 



Figs. 1-6 show variation of t a h — to of eq.(13) for the variation of fco(= go/Do) and ki(— gi/Di). Figs.l, 3 and 5 
correspond to D\ > for n — 4, 12 and 25 respectively while Figs. 2, 4 and 6 correspond to D\ < for n — 4, 12 and 25 

respectively. 



(/(r) = 0) with I = (n + l)/2 will be considered (as in the previous section). Now choosing the initial time 
ti = 0, the explicit solution for i?(t, r) can be written as 



t(r) = 



V2n 
30) 



D{r) 



n + 1 



5 (r)i?(™+ 1 )/ 2 - \ D(r) 



n+1 



g ( r ) r (™+ 1 )/ 2 



(15) 



Thus the expression for the singularity time for the shell of radius R is given by (R(t s (r), r) = 0) 



ts(r) 



9{r) 



JDM- \ D(r) — g ( r ) r (n+i)/2 

V n + 1 



(16) 



and consequently the time for central singularity is 



to 



2/i 



'D n -i D n - 



3" 



go y v n + 1 

Here the polynomial form of D(r) and g(r) are taken in the form 

D(r) = D Q r n+1 + D k r k+n+1 
9{r) = 9o r (n+1 V 2 + g . W+(" +1 )/ 2 



(17) 



(18) 



where D ,3o are constants and D k (< 0) and g .(< 0) are the first non-vanishing term beyond Do an d go 
respectively. Now using these expression for D(r) and g(r) the time of collapse of a typical shell of radius R 
(i.e., t s (r)) becomes (see eq (16)) 



f^D k I 1 
t s {r) =t +\ - — 



2 3o \ VAi 



Do - ^30 . 



3o 



^+1)^0-^30 



'A. - A /I>o-^flD 

3o 



(19) 

Now the equation of the outgoing radial null geodesic (ORNG) which passes through the central singularity 
can be chosen to be (near r — 0) 



toRNG = h + a r ? 
where a(> 0) and £(> 0) are constants. 

In the polynomial form for D(r) and g(r) in the equation (18) one can choose two possibilities: 

(i) k <j , (ii) k> j 

Case I : k < j : 
Here near r = 0, the expression for t s (r) can be written as (see eq.(19)) 



(20) 



t s (r) =t a - 



n Dk 

2 3o 



1 



Do 2 



n+l 



</n 



iDo 



r k , (D k < 0) 



(21) 



In order that null geodesic passes through the shell of radius R before the trapped surface is formed there, 
comparing (20) and (21) one gets (toRNG < ts(r)) 



Dk rn 



(a) £ > k or (b) £, = k and a < 

go V 2 

When £, > k then near r = the solution for R simplifies to 



1 



1 



A, 



n+l 



30 



^Do 



(22) 



R = r 



n+l 
An 



got 1 



n+l 



2n 



t A D 



n+ 1 



3o 



2,/Dn- 



71+ 1 



3o 



2/n+l 



(23) 



Further for the given metric an ORNG should satisfy 



dt 

dr 



R' + Rv' 



(24) 



Now using (20) and (23) in (24) one gets (up to leading order in r) 



a^- 1 = 1 



(n + l)D k to 



2V2^JD Q - ^go 



This gives 



e = i 



2k 

- 1 + ^TT 



2k 1 / 2fc 

> and a — — 1 + v . -\ 

n+1 £ V n+1 



2/(n+l) 



r— , (i/_ t ? 0) 



(25) 



(n + l)D k t 



2V2n J D ^.go 



2/(n+l) 



(26) 



As £ > fc, so from the above relations (26) 



k < 2+1 and £ < 2±f, 

n— 1 ^ n— 1 ' 

i.e. one could have (n > 2) 
fc = l, £ = s±? 
On the other hand for £ = fe, as before A; = ^3^ an d 



(27) 



n — 1 

n+1 



n+1 20^0 / 2 

4 I n V n V n+1 



,9o 



(1 



2ag t 



2a J- 
V n 



Ai- 



-.90 



A* yi y 


(l-n)/ 


(n+1) 




f A*o | 


^0- 


^+T50 


jj 



(28) 



Thus for £ > k, k has only one value (namely k = 1) and n can take any value (> 2) while for £ = fc 
the only possible values of n are 2 and 3 only. Therefore, geodesic equations are possible in any dimen- 
sion for £ > k but it is only possible up to five dimension for £ = k. In other words, for £ = k, naked 
singularity is possible only up to five dimension which supports the results in the previous section. Lastly, 
it should be mentioned that the other choice namely k > j is similar to the above and hence not presented here . 



IV. DISCUSSIONS AND CONCLUDING REMARKS 



A detailed analysis of the central curvature singularity as the final state of collapse in the (n + 2)-dimcnsional 
quasi-spherical Szekeres' model has been done for matter with anisotropic pressure (i.e., both radial and 
tangential pressures are non-zero and distinct). The local visibility of the central singularity has been discussed 
by comparing the time of formation of trapped surface and the time of formation of central shell focusing 
singularity while global visibility is examined by considering only the radial null geodesies (for simplicity). 
Most of the results are very similar to that for four dimension in ref. [13]. It is to be noted that though we 
have considered quasi-spherical Szekeres model but still these are valid for TBL model. In fact, in Szekeres' 
solution if we assume the function v to be independent of r (i.e., v' = 0) then Szekeres' model can be converted 



to TBL model by the following coordinate transformation: 



x\ = Sin$ n Sin$ n -\ SinQ^Cot^Bi 

X2 = Cos9 n Sin8 n -i SinQ^Cot^Bi 

£3 = Cos$ n -\Sin6 n -2 SinQiCotiflx 



x n = CosQiCot\Qx 



Further, the radial pressure is assumed to be a function of 'r' and '£' only with the form p r = g(r)/R l . 
Throughout the paper I is chosen to be less than n + 1 (i.e., I < n + 1) [It is to be noted that for I = n + 1 
due to appearance of a logarithmic term the calculations become very much complicated while for I > n + 1 
the results are not of much interest]. It has been shown that for 31 = = D± i.e., for hi = 0, naked singularity 
is only possible up to five dimension — a result identical in dust collapse. The above choice (i.e., hi = 0) gives 
po = (n+ l)/io for i/_i > (—1) and pi = for i/_i > — 1. Hence for v_i > — 1, if the initial density gradient falls 
off at the centre (r = 0) then one can say that six dimension plays as critical dimension for naked singularity, 
a distinct result in higher dimension. Further with the choice I = (n + l)/2, a detailed comparative study 
has been done between the time of formation of trapped surface and that of central singularity and Table I 
shows all possibilities for the parameters involved in the expression. Also in the figures 1-6 we have shown 
graphically the time difference t a h — to for 6, 14 and 27 dimensions for Di > and Di < respectively. 

As in dust collapse, in this case we have definitely a black hole (or naked singularity) if the initial density 
gradient at the centre is positive definite (or negative definite). But in the indefinitcness in the sign of pi 
one may note that if the initial density and radial pressure has identical behaviour (i.e., increase or decrease 
simultaneously) then even with initial negative density gradient (at the centre) it is possible to have black hole 
as the end state of collapse, while if the initial density and pressure have opposite nature (i.e., one increase when 
other decreases and vice versa) then the behaviour is identical to dust collapse. Therefore one may conclude 
that pressure tries to resist the formation of naked singularity. 

APPENDIX I 

Initial hypersurface and the physical parameters: 

Suppose the collapsing process starts on the initial hypersurface (t — ti) and we have R = r there. Then the 
expressions for energy density, radial pressure and the tangential pressure, at the beginning of the collapse are 



/ \ /. \ h' (r) + (n+l)h(r)v f 

Pi\'i *^1 7 ■•■7 •^n) PV'ii'^i'^li ■•■; %n ) r"(l+ri/') 

(20) 
p (r,xi,...,x n ) =p T (t = U) = Sip. 



n{l-\-ri/ / ) 



9'(r) 



nr l 1 (1+ri/ 7 ) 



where h(r) = H{r,U) = D(r) - ^L. 



9(r) r n-l+l 



For smooth initial data h(r) and g(r) to be C°° functions and hence one can choose the following series 
expansions 



p l {r, x x , ...,x n ) = Ejlo Pi r ° ( 30 ) 



v'{r,x 1 ,...,x n ) =Y!jL-i v j r3 > i v -i > - 1 )- 

p (r,xi,...,x„) =E^oft' ^ 



It is to be noted that in the above series expansions the coefficients Dj's and gj's (j = 0, 1, ...) are purely 
constants while Pj's, i/j's and Pj's are functions of ar^'s (i = 1, ..., n). Also these coefficients are not independent 
but are related among themselves through the relations in eq (29) as follows: 



For v , > — 1: 



P0 = <?0, Pl=.9l{l+ n(1+ 1 „_ i) }, P2 = .92 {l + XT+^) } - n(i+,T> ' - 
/ L1 U (n+2) + (n+l)7/_ 1 ( n+ 3)+( n+ l) ! ,_ 

po = (n + l)/io, Pi = 1+7^ &i> P2 = i+t^ ft 2 - (!+„_[)* , 



(31) 



For i/_j = — 1: 

Po = 9o + ^, Pi=5i(l-^%)+^, P 2 =. 92 (l-2^)+M_|-) 5l + |^, 

Po = (n+l)/*> + £, Pi = ^f + ^i((n + l)-^), p 2 = ^ + ft 2 ((„ + l)-^)+ft 1 I^5l, . 
with hi = Do %7T, i = 0, 1, 2, .... 

APPENDIX II 

Solution of the evolution equation (11) with f(r) = 0: 
For the initial choice R = r at t = t%, the explicit solution is 



(32) 



V2^r— 1 ff (r)r™ +1 -' V2^ i?— 1 ff(r-)i?" 

^^"(n + l)/^ 2Fl[ 2' 6 ' 5+1 '^)(n + l-0 F (n + l)/DW^ l[ 2' 5 ' 6+1 'i?(r)(n + l-0 ] (M) 

where 6 = 2 "^ and 2-Fi is the usual hypergeometric function and I ^ n + 1. 

If i = t s (r) stands for time of collapse of the r-th shell i.e., R(t s (r), r) = then we have 



i/2r7 r(" +1 )/ 2 1 rjfrV n+1- ' 

t.(r)-ti= / lH " -— ^[1,6,6+1 g \7 ] (34) 

(nH-l)-vAD(r) 2 D(r)(n + 1-Z) 

Note that £ = t s (r) is a monotonic increasing function of r (i.e., t' s (r) > 0) for the shell focusing singularity. 
Using equation (12) 

fa,l(rM ^(n + l)V^M 2 l[ 2'' + 'W(n + l-0 ] (n + l )v ^M ^ [ 2' M+1 ' Z?(r)(n + 1-1) 

(35) 



10 

A comparison of equations (34) and (35) shows that the shell focusing singularity that appears at r > is in 
the future of the apparent horizon. However, the time of occurrence of central shell focusing singularity (which 
is of main interest here) is given by 



to = lira t s (r) 
r^Q 



^(^Sk^liM + M, (* 



(36) 
go 



D (n+l-l) 



where the series form of g(r) and D(r) (from cq. (30)) have been used in evaluating the limit. Thus for the 
restriction I < n + 1 , the explicit form of the difference between t a h(r) and £q is 



t a h(r) -t = 



Fh~ 17 D 3/2 D 1 1 (gpgi - gigo) 5/2 3 ] 

V2[\'^Tl~ 2fl[ 2' b ' b+1 ' z]+ (n+l-0(3 + 3n-20 ^ 2^-, 6+ 1,6 + 2,,]^ 



>+i 



2^n^ J D - 1 PI l llM /2Do 

\ n 



0(r») - Tr ^ ^[-,6,6+1 



z 



(37) 
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